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ABSTRACT 

One of the main endeavors of fundamental astrometry is to establish a prac- 
tical realization of a non-rotating, inertial reference frame anchored to celestial 
objects whose positions are defined in the barycentric coordinates of the solar 
system matching the current level of astrometric observational accuracy. The de- 
velopment of astrometric facilities operating from space at a microarcsecond level 
of precision makes the non- uniformity of the galactic motion of the barycenter 
an observable and non-negligible effect that violates the desired inertiality of the 
barycentric frame of the solar system. Most of the observable effect is caused by 
the nearly-constant (secular) acceleration of the barycenter with respect to the 
center of the Galaxy. The acceleration results in a pattern of secular aberration 
which is observable astrometrically as a systematic vector field of the apparent 
proper motions of distant quasars. 

We employ the classic approximations of planar epicycle and vertical har- 
monic oscillation for the Sun's galactic motion to estimate the magnitude of 
secular acceleration components in the galactic coordinates and show that these 
approximations are adequate for the SIM space mission. We employ the vector 
spherical harmonic formalism to describe the predicted field of proper motions 
and evaluate the amplitude of this field at each point on the celestial sphere. 
It is shown that the pattern of secular aberration is fully represented by three 
low-order electric-type vector harmonics, and hence, it is easily distinguishable 
from the residual rotations of the reference frame and other possible effects, such 
as the hypothetical long-period gravitational waves, which are described by other 
types of vector or tensor harmonics. Comprehensive numerical simulations of the 
grid astrometry with SIM PlanetQuest are conducted assuming that 110 optically 
bright quasars are included as grid objects and observed on the same schedule 
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as regular grid stars. The full covariance matrix of the simulated grid solution 
is used to evaluate the covariances of the three electric harmonic coefficients, 
representing the secular aberration pattern of proper motions. We conclude that 
the grid astrometry with SIM PlanetQuest will be sensitive to the main galac- 
tocentric component of secular acceleration, arising from the circular motion of 
the Local Standard of Rest (LSR) around the galactic center, while the peculiar 
acceleration of the Sun with respect to LSR is expected to be too small to be 
detected with this astrometric space interferometer. 

Subject headings: astrometry — techniques: interferometric — reference systems 
— Galaxy: kinematics and dynamics — quasars: general — dark matter 

1. Introduction 

The Sun accelerates as it travels in the galaxy. The acceleration is small (approximately, 
6 mm s~^ yr~^), and it takes millions of years to change the Sun's velocity vector significantly. 

The instantaneous velocity of the Sun in a non-rotating galactic reference frame with the 
origin at the centcr-of-mass of the Milky Way causes all observed directions to both galactic 
and cxtragalactic sources to deflect from their "true" directions in a systematic way. This 
phenomenon has the same spccial-rclativistic nature as the annual stellar aberration induced 
by the orbital motion of the Earth around the Sun. Due to the aberration the observed 
position of a fight source is displaced toward the direction of the instantaneous velocity of 
the observer with respect to an inertial reference frame at rest by an amount proportional 
to the velocity magnitude (Kovalevsky & Seidelmann 2004). 

For an observer located at the barycenter of the solar system, the instantaneous effect 
of the relativistic aberration due to the galactic motion of the solar system is not directly 
observable at any fixed instant of time because the velocity-induced aberration pattern is 
constant. But since the motion of the solar barycenter with respect to a fixed, non-rotating 
galactic reference frame is non-rectilinear, the quasi-stationary (secular) change of its velocity 
vector causes the aberration pattern to change on the sky gradually as time progresses. 
Given astrometric measurements at a microarcsecond (yuas) level of accuracy conducted over 
a time span of at least several years, the changing aberration pattern can be observed in 
two ways. An astrometric facility, such as Very Long Baseline Interferometry, capable of 
measuring long arcs between light sources at two separate epochs, could directly detect and 
measure the distortion of the grid of a global astrometric frame based on quasars taken as 
reference objects (Fomalont & Reid 2004). This distortion would be a direct consequence 
of the directional difference in Sun's galactic velocity between the two epochs, as well as 
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other possible phenomena, e.g., unconstrained residual rotations of the reference frame, 
gravitational lensing from invisible dark matter (Liddle 1999), long-periodic gravitational 
waves of cosmological origin (Pyne et al. 1996; Gwinn et al. 1997), etc. 

A more tenable approach to detecting the secular aberration is to observe the global 
vector field of the proper motions of quasars provoked by the slowly changing aberration 
pattern that is induced by the secular acceleration of the Sun (see Fig. 1). This idea has 
circulated among astromctric community for a number of years (Ferryman 2000). Nowadays 
it becomes technologically feasible and can be rendered within the framework of the existing 
projects of space astrometry which are currently under development in the national (SIM 
space mission) and European space agencies (Gala space mission). 

Froper motions of stars projected on the celestial sphere are customary expressed in 
angular units (e.g. arcseconds) per year, and stem from the spatial motion of stars with 
respect to the solar barycenter. Fhysical proper motion of each star is proportional to the 
tangential, that is orthogonal to the line of sight, component of the spatial velocity and 
inversely proportional to the distance to the star. For this reason, physical proper motions 
of quasars, which are separated from the Sun by distances of many Megaparsecs, are orders 
of magnitude smaller than 1 //as yr~^, and, thus, will be hardly observable in the near future. 
The secular aberration, however, mimics the physical proper motion of a distant celestial 
object and makes all of them, including "infinitely" distant quasars, move on the sky with 
observable proper motions of ^ A yuas yr~^, which can be detected with current technology. 
Measuring the secular aberration is vitally important for evaluating the total mass of the 
Milky Way including both visible and invisible matter in its bulge and halo which will provide 
us with a precise estimate of the amount of dark matter in our Galaxy. 

This paper discusses the astrometric effects of secular aberration which can be detected 
by the SIM space mission. The basic relationship between the galactocentric acceleration of 
the solar system and the observed pattern of proper motions of distant quasars is established 
in section 2. The planar epicycle theory and a harmonic vertical-motion approximation are 
used in section 3 to evaluate the peculiar acceleration of the Sun with respect to the Local 
Standard of Rest (LSR). A decomposition of the resulting proper motion field in orthogonal 
vector spherical harmonics is given in section 4. The vector-harmonics technique and the 
simulated global covariance matrix of a number of quasars potentially observed with the 
Space Interferometry Mission are used to compute the expected accuracy of measuring the 
secular acceleration for this project (section 5). 
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2. Proper motions and seculeir aberration 

Let us choose baryccntric Cartesian coordinates (X, F, Z) such that the plane of the 
galaxy is the X — Y plane and the X axis is directed towards the center of the galaxy. 
The axis Y is aligned with the direction of the galactic rotation. Let us introduce spherical 
coordinates (r, /, b) such that the unit basis vectors of the two coordinate systems are related 



by the following equations (Binney & Merrifield 1998) 

Cr — cosfccosZ ex + cos6sinZeV + sin6ez (1) 

ei = — sin / ex + cos/ ey (2) 

Cb = — sin b cos I ex — sin b sin I eV + cos b ez , (3) 



where the angular coordinates (/, b) are the galactic longitude and latitude respectively, and 
the direction to the galactic center is given by the unit vector ex- The geometric direction 
to a star at coordinates (r, /, b) is given by the unit vector e^ = K. Let us assume that at a 

— * 

given epoch to the orbital velocity of the Sun around the galactic center isV— {Vx, Vy, Vz) 
and the acceleration is ^4 = {Ax, Ay, Az). 

Observed direction to the star measured by a moving observer differs from its true 
geometric position, because of the stellar aberration that is proportional at each instant of 
time to the velocity of the observer with respect to the static frame. Since the velocity vector 
of the solar system is not constant due to the galactocentric acceleration, the aberration angle 
changes progressively. Therefore, at another epoch t, the star is seen by a fictitious observer 

— * 

located at the solar barycenter in the direction, k, given by 

k^K^^Kx{V xK)^-^Kx{AxK){t-tQ) , (4) 

where " x " denotes the Euclidean cross-product of two 3-vectors, c is the speed of light 
in vacuum, to is the initial epoch, t is the time of observation, and we assume that the 
contribution of the first and higher order derivatives of acceleration is neghgible. The first 
two terms in the right-hand side of Eq. (4) are constant if one neglects the secular parallax 
(Binney & Merrifield 1998), which is irrelevant for quasars and will not be considered in 
this paper. The third term in the right-hand side of Eq. (4) yields the secular aberration 
that makes all objects beyond the boundaries of the solar system to move in the sky with a 
proper motion 

fl = IJ,i ei + fib Cb , (5) 

where the longitudinal and latitudinal components of jl are proportional to the acceleration 
of the solar system projected on the celestial plane 

III — - {—Ax sin I + Ay cos I) , (6) 
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fj,b — - {—Ax sinbcosl — Ay smbsml + Azcosb) . (7) 
c 

If one assumes that the acceleration of the solar system barycenter has only Ax = A compo- 
nent, the equations (5)-(7) for the effective proper motion caused by the secular aberration 
are simplified and reduced to 

jl— (sin / + sin 6 cos / Cb) . (8) 

c 

Proper motion vectors /2 of a given number of objects represent a discrete vector field on the 
sphere which can be decomposed in a set of vector spherical harmonics (Thorne 1980). The 
largest galactocentric component of the secular aberration can be determined from global 
astrometric observations as a systematic dipole component of this vector field. For quasars 
the problem of its determination is simpler than for stars since they have negligibly small 
proper motions caused by their peculiar velocities with respect to the Hubble flow. Therefore, 
the secular aberration can be directly measured from the observed proper motions of quasars. 

The magnitude of the secular aberration effect in the case given by Eq. (8) is (Ferryman 
2000) 

= - sin C , (9) 
c 

where sin ( — y/1 — cos^ I cos^ b and ( is the angle between the direction towards the galactic 
center and that to the star. In what follows, we investigate a more accurate approximation 
of Eq. (8) that includes all three components of the Sun's acceleration, and evaluate the 
effect of the secular aberration more adequately in terms of vector harmonics. 



3. Circular and peculiar accelerations of the Sun 

The velocity vector of the Sun in the galaxy is commonly considered to be the sum of two 
components (Binney & Merrifield 1998). The first (and largest) component is the motion 
of the so-called. Local Standard of Rest (LSR). By definition (Binney & Merrifield 1998), 
the LSR is involved in a circular planar motion around the center of mass of the galaxy at a 
constant rate with a period Pq. The second component of the solar velocity is the differential 
(peculiar) motion of the Sun with respect to the LSR, which makes the Sun's orbit to be 
non-circular and non-planar. We adopt a solar peculiar velocity of Vq = (10.0,5.3,7.2) 
kms~^ as given by (Dehnen & Binney 1998). The absolute value of the rotational velocity 
of the LSR is 

(^) (^) 
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where Ro is the distance of the Sun from the galactic center and the basic quantities Pq and 
i?o are normahzed to their best known values (Binney & Merrifield 1998). 

The galactocentric acceleration resulting from this circular motion is 

and the maximum value of a proper motion caused by this acceleration is 

As follows from equation (9) the maximum proper motion is achieved for objects lying on 
the great circle orthogonal to the direction to the galactic center. The proper motion vectors 
from this largest component of the LSR acceleration are directed toward the galactic center 
(see Fig. 1). 

In order to estimate the magnitude of the secular aberration caused by the peculiar 
acceleration of the solar system with respect to the LSR, we employ the classic epicycle 
approximation for the planar motion of the Sun, and a harmonic force approximation for its 
vertical oscillation about the galactic plane (Marochnik & Suchkov 1984; Binney & Merrifield 
1998). The peculiar velocity components in this model are (Makarov et al. 2004): 

kA ■ ■ k 

X — — Xq— sin kt + Xq cos kt + Yq-— sin kt , (13) 
B 2B 

2A(A — B) k 1 

Y = -Xo — ^— '-{l-coskt)-Xo— sin kt + Yo— {A- {A- B) cos kt) , (14) 

B 2B B 

Z = —ZQUsinut + Zocosi't, (15) 



where the planar epicycle frequency k — ^/—4:B(A — B) is expressed in terms of the Oort's 
constants A and B, [Xq, Yq, Zq) are the galactic coordinates of the Sun with respect to the 
LSR at time to, and {Xq^Yq^Zq) are their velocity components. The vertical frequency, 
V = 2Ti/Py, of the oscillatory motion of the solar system about the galactic plane is an 
ill-constrained parameter that depends on the vertical period P^, which we assume to be 
60 Myr. The harmonic vertical oscillation is a linearized approximation, which ignores the 
second-order terms in the expansion of the vertical component of the gravitational force over 
coordinate Z. According to the current knowledge of the mass distribution in the Milky Way, 
this second-order terms are less than 3% of the main harmonic term (which is proportional 
to Z) within 70 pc of the plane (Makarov et al. 2004). Hence, they can be neglected in our 
estimations. The harmonic acceleration derived from the equations (13)-(15) is accurate to 
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about 0.05% for the Sun, which is lying at 9 ± 2 pc above the galactic plane (Marsakov & 
Shevelev 1995). 

The epicycle theory is also based on a linear approximation derived from the expansion 
of the local horizontal component of the gravitational force over radial distance from the 
LSR R — Rq. The non-linear terms in the expansion of the horizontal force can reach 15% 
of the magnitude of the linear term at a distances of 1 kpc from the LSR. The Sun does not 
stray farther than about 300 pc from the LSR in this direction. Our conjecture is that the 
quadratic terms in the radial-force expansion will be less than 1.5% of the epicyclic (linear) 
term. Because the quadratic terms correspond to the time derivative of the solar peculiar 
acceleration, they appear to be negligibly small. Little (if anything) is known about the 
non-radial local gravitational forces resulting from the influence of the galactic bar, spiral 
arms and other possible massive attractors including dark matter. We do not attempt to 
include them in the calculations given in the present paper but a theoretical study of their 
impact on the peculiar components of the solar system motion would be highly desirable. 

Differentiating equations (13)-(15) and setting t — and Xq — yields the solar pecuhar 
acceleration 

ax = -2{A-B)Yo, (16) 
ay = 2{A-B)Xo, (17) 
az = -Zou^ . (18) 

The corresponding numerical estimates of the peculiar acceleration components are, respec- 
tively 

ax = -0.27-10-6 ( ^^i— ^ -] ( [kms-' yr-^l (19) 

V0.025 kms-^pc-V \^5.3kms-^y ^ ^ ^ ^ 

ay - 0.51 ■ 10-6 ( ^ ( ^ A [kms-i yr-^l (20) 

V 0.025 km s-^pc-V ^lO.Okms-^y ^ ^ ^ ^ 

Equations for the maximum proper motions corresponding to this peculiar acceleration ex- 
pressed in //asyr"^, are obtained by dividing equations (19)-(21) by 1.45 • 10"^^ which is 
a conversion factor of units used in the calculation. We find that the residual proper mo- 
tion field caused by the peculiar acceleration of the solar system with respect to the LSR 
is smaller in amplitude than 1 /xasyr"^, that is about 10 times smaller than the secular 
aberration (12) caused by the galactocentric acceleration of the LSR shown in equation (11). 
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4. Vector harmonic analysis of stellar proper motions 



A global discrete pattern of proper motions is a vector field /!(/, b) on the celestial sphere 
which can be expanded in orthogonal vector functions of spherical coordinates 



oo 



(22) 



where Tj and Sj are orthogonal vector harmonics, called magnetic and electric harmonics 
respectively (Thorne 1980), and tj and Sj are the coefficients of the expansion. The vector 
harmonics Tj and Sj can be expressed in terms of the partial derivatives of the scalar spher- 
ical harmonics Vj (Arfken & Weber 1995) with respect to galactic longitude and latitude. 
Specifically, one has (Vityazev & Shuksto 2004) 



dVAl,b) 
db 



1 dVj{l,b) 
cos b dl 



eb 



Sj{l,b) = Q 



1 dVj{l,b) ^ dVj{l,b) ^ 

~l~ ryj ^6 



cos b dl 



db 



(23) 
(24) 



where Cj = Cnm are constant normalization coefficients making the harmonics orthonormal, 
and the cumulative index j = {nm} counts the spherical harmonics Vj = Vnm over orders 
n = 0, 1, . . . , cxD and degrees m = 0, 1, . . . , n. We do not specify the normalization constants 
Cj because real observations will provide a vector field of proper motions sampled at a 
number of discrete points on the celestial sphere corresponding to the observed quasars. 
This constant is determined in the process of calculation of the covariance matrix of vector 
harmonic coefficients, as explained in the following text. 

The scalar spherical harmonics are given in terms of the associate Legendre polynomials 
as follows (Arfken & Weber 1995) 



Km = \^ , , , Pr(sm6)cosm/, 



271 (n + rn)V 
'2n + 1 (n -m)! 



P^{smb) sin ml . 



(25) 
(26) 



27r {n + m)\ 

A useful identity for differentiating the Legendre polynomials is 

^ = ^[{n + 1) sin6P™(sin6) - (n - m + 1)P™ i(sin6)] , (27) 
and it will be used in order to derive the vector spherical harmonics in equations (23)-(24). 
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The first set of spherical harmonics generating a pair of vector harmonics, are the zonal 
harmonic Vio ~ sin 6 and the sectorial harmonics V^^ ~ cos 6 cos Z, V^^ ~ sin 6 sin L The 
vector harmonics corresponding to Viq are 

Tio(/,&) = Ciocosfeel (28) 
Sio{l,b) = Cwcosbel. (29) 

The vector harmonics corresponding to Vfi are 

T^i = Cu (cos / sin 6 6/ — sin / e^) , (30) 

— * 

— Cii(cos/sin6 Cfe + sin/ e;) , (31) 
and those corresponding to V^i are 

= Cn (sin / sin 6 q + cos / el) , (32) 

— * 

S^-^ — Cii(— cos/ el + sin/sin6 Cb) . (33) 

Numerical factors Cio and Cn in Eqs. (28)-(32) are constant normalization coefficients 
which are not required in our discussion. 

The vector harmonic decomposition of the proper motion field defined by Eq. (5) can 
be given now in the following form 

fla^l {-AxS'n - AySt, + AzSto) , (34) 

which chicidates that the proper motion field caused by the secular aberration is represented 
by the three low-order electric dipole harmonics of first order (n = 1) and it is independent 
of the magnetic-type harmonics. This systematic vector field affects all observable directions 
and can be extracted from the randomly distributed vector field of stellar proper motions. 
However, this subtle effect is submerged in the pattern of the physical proper motion field of 
galactic stars, which is roughly three orders of magnitude larger than the secular aberration 
effect (Oiling & Dehnen 2003). For quasars, because of their negligible physical proper 
motions, the systematic pattern of secular aberration stands out clearly, and is a dominant 
component of the quasar proper motions. Fig. 1 shows the sky distribution of the proper 
motion pattern caused by the acceleration of the solar system barycenter. 
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5. Numerical Simulations of the Secular Aberration and Accuracy of Its 

Measurement by SIM 

5.1. SIM Astrometric Grid Model 

SIM PlanetQuest is a NASA Origins project dedicated mostly to the search of planets 
around nearby stars by astrometric interferometry As a necessary condition of achieving this 
goal, SIM will create a global astrometric reference frame over all the sky to an unprecedented 
accuracy of 4-5 //as in position, based on wide-angle measurements of about 1300 reference 
grid stars. It is crucial for the astrometric coordinate grid to include a sufficient number of 
optically bright quasars, which constrain the parallax solution of the astrophysical targets 
in such a way that the large-scale random distortions of the parallax error distribution are 
drastically reduced (Makarov & Milman 2005). The same grid quasars can be directly used 
to determine the secular acceleration of the Sun from measuring the secular aberration effect 
as explained in the previous sections of this paper. 

To estimate the measurement accuracy of the secular aberration effect for SIM Plan- 
etQuest we make use of full-scale numerical simulations of SIM observations and the astro- 
metric grid reduction model (Makarov & Milman 2005). The simulation model is rather 
sophisticated and realistic in that it includes the SIM instrument model and incorporates 

the specific technology of SIM astrometric measurements. The grid reduction model has 
about 160,000 free fitting parameters (unknowns) and deals with nearly 315,000 numerically 
simulated delay measurements taken during supposed 5-year mission life span. The adopted 
one-step solution of the global astrometric grid is found by the least-square method which 
produces the covariance matrix of fitting parameters. Among the original 160,000 parame- 
ters only 6,500 determine the astrometric grid, while the rest are field- dependent instrument 
and baseline-orientation parameters ^. In this paper, we use only that part of the covari- 
ance matrix which contains astrometric information about the grid reference quasars directly 
related to determination of the secular aberration effect. 

Let us define a column array of vectors = {/Ii, /I2, • • • , ySw} comprised of the proper mo- 
tion vectors of N quasars. We define a matrix S representing A x 3 array, whose columns are 
the vector harmonics taken at the points 1, 2, . . . , A. Let us also introduce a column vector 
a = {ui} consisting of three coefficients ai, 02 and 03 representing respectively Ax, Ay, Az 
components of the secular acceleration vector (see equation (34)) which is to be determined. 
These parameters are found in our model by the least squares method employed for solving 



^See (Makarov & Milman 2005) and references therein. 
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the system of linear equations 

S • a = /u , (35) 

where the dot "•" denotes the inner product with respect to the index i — 1,2, 3. If a is a 
least squares solution of equation (35), its covariance is (Bard 1974) 

Cov[a] = E [dd^] = {S^Sy' (S^Cov[iu]S) (S^S)"' ; , (36) 

where E[r] denotes the expectation value of random variable r, S-^ is the transposed matrix 
S, and Cov[/Li] = E[/li • /x^]. We emphasize two points: 

1. Solution d gives the components of the secular acceleration expressed in the same 
dimensional units as the observed proper motion vector field /i, that is in ^as yr~^. This 
choice of units determines coefficients Cj in equations (23)-(24) which are naturally 
fixed for a given discrete sample of points by the normalization of the matrix S in 
equation (36). 

2. The vector spherical harmonics used for computing the fitting parameters are not 
orthogonal because they are sampled at a finite number of points on the celestial sphere. 
This accounts for the presence of terms S^S in equation (36). 

5.2. CovEiriance Analysis 

Covariance Cov[/Lt] is a x matrix which is easily calculated from the full covariance 
matrix of the astrometric grid solution (Makarov & Milman 2005) in the ecliptic coordinate 
system. In doing this, we do not have to rotate the covariance matrices from the ecliptic to 
the galactic coordinates because the inner products (/Xfc-/Xm), for which the covariance matrix 
in equation (36) is computed, are invariant with respect to rigid rotations of the coordinate 
system. However, the proper motion vectors should be projected onto the galactic coordinate 
axes by making use of an appropriate orthogonal transformation (e.g., ESA 1997, Vol. 1) 
before solving the system (35), because the vector harmonic functions are not invariant with 
respect to rotations of spherical coordinates. We also notice that for any pair (k, m) of the 
grid quasars the expected value of the inner product of their proper motions, Covf/Li]^^, is 
the sum of the two corresponding covariances of the coordinate proper motion components. 

Following this algorithm, we carry out a full-scale numerical simulation of the secu- 
lar aberration problem, using the grid solution method described in detail in (Makarov & 
Milman 2005). We assumed in this paper that a set oi N — 110 quasars having a fairly 
uniform distribution on the sky, is included in the SIM reference grid. The astrometric frame 
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consists of a sufficiently dense grid of reference stars and is required to carry out astrometric 
reductions of the science target data collected in the wide angle regime. The need for a 
reference grid arises from the fact that the space interferometer will measure the optical 
path delay of the light-wave front from celestial sources with respect to its own baseline 
vector, which has to be determined from the same observations in a self-consistent manner 
^. The majority of the SIM grid objects will be distant red giant stars of rriv — 11 mag- 
nitude, spectroscopically vetted as non-binary stars. However, it is imperative to include 
in the grid a number of optically bright quasars. Having parallaxes much smaller than 1 
lias, the quasars will serve to constrain the global parallax solution significantly reducing 
the zero-point and other large-scale parallax errors for scientific targets. Positions of the 
grid quasars determined with SIM will be used to remove the residual distortions of the grid 
reference system based on the red giants, with respect to the extra-galactic reference system 
based on the quasars (Johnston et al. 2003). We note here, that in order to make the SIM 
reference frame truly inertial to accuracy better than 1 yuas yr~^, the residual spin present 
in the proper motion field and caused by the indefinite rotation of the grid solution, must 
be removed as well. Rigid rotations (spins) of the reference system around the coordinate 
axes X, Y and Z are represented by the first three magnetic- type vector harmonics of the 
observed proper motion field (Vityazev & Shuksto 2004). Hence, residual spins are orthog- 
onal to the secular aberration effect and, for this reason, can be neatly separated from it 
in the vector harmonic space. This point can be effectively demonstrated by comparing the 
pattern of the proper motion field induced by the secular aberration and shown in Figure 
1, against the patterns of the proper motion fields induced by residual spins of the global 
solution about X, Y, Z axes shown in Figures 2, 3, and 4 respectively. 

The covariance matrix. Gov [a], of the secular acceleration components is given in Table 1. 
The diagonal elements of the matrix are the unit-weight variances of the corresponding 
X, Y, Z components of the proper motion field (34). For example, the standard deviation of 
ai in units of //as yr^^is estimated as follows 

a[ai] = VCovii[a](7o = 0.076 (7o, (37) 

and similarly for other coefficients, where is the standard error of a single measurement 
(the interferometer path delay), expressed in /xas. For quasars as bright as niy = 14 magni- 
tude, single measurement errors of 15-i-20 //.as could probably be achieved without spending 
too much integration time on these objects. Hence, the largest component of the proper mo- 

— * 

tion field in S^-^, which is proportional to the galactocentric component Ax of the solar 



^Similar problem exists in radio intaferometry (Sovers, Fanselow & Jacobs 1998). Its solution also relies 
upon self-consistent construction of a fundamental reference frame. 
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acceleration, will be determined with an error, a{Ax) — 1.5 //as yr""^as follows from equa- 
tion (37). This accuracy is sufficient for a statistically robust detection of the acceleration 
at a signal-to-noise ratio of roughly 3. 

6. Discussion 

We have seen that, upon observing 110 optically bright quasars as grid objects, SIM 
PlanetQuest will be able to detect the major galactocentric component of the acceleration 
of the Sun. At the same time, the anticipated accuracy of SIM is not high enough to detect 
the peculiar acceleration of the Sun with respect to the LSR, nor to refine the fundamental 
constants Pq and Rq. In order to do this a dedicated interferometric mission is required. 

Proper motions of the grid stars and the science target stars will also include the sec- 
ular acceleration components. These effects will be swamped by the stronger patterns of 
differential galactic rotation, asymmetric drift, local and large-scale dilation, disk warp, per- 
turbations from the spiral arms and large streams of galactic stars, so that stellar proper 
motions will probably be useless to study the subtle secular aberration effects. The latter, 
nonetheless, have to be removed from the observed proper motions of stars to achieve the 
highest accuracy and unbiased estimates in the galactic dynamics studies. By the same 
argument, the residual spins of the stellar proper motion field (which, in fact, may be rather 
large) can be reduced to ~ 1 //as yr~^ accuracy from quasar measurements. 

Rigid rotations (spins) of SIM's reference frame are parts of the null space of a global 
astrometric solution (Lattanzi, BucciareUi & Bernacca 1990; Makarov & Milman 2005), 
and, therefore, can not be determined from observations with SIM. They are confined to 
three low-order magnetic vector harmonics, while the vector field of proper motions caused 
by the secular acceleration is represented by the three electric harmonics of zero and first 
order. SIM will establish an optical reference frame based on a number of reference stars. 
In order to make this frame as inertial as possible observation of bright quasars are of 
utmost importance for making the residual rotations of this frame appropriately small. It 
is important to ensure that the grid quasars are uniformly distributed on the sky and each 
of them is observed to roughly the same precision, to keep the correlations between the 
astrometric parameters as small as those calculated in our numerical simulations (Table 1). 
These requirements prevent the random errors from propagating between the constituents 
of the vector harmonic space. Higher-order harmonics of the quasar proper motion field will 
also be of a vast interest to examine since they will provide an independent estimation of 
the external accuracy of SIM. 
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In principle, some additional large-scale pattern in the observed proper motions of grid 
quasars may be induced by the non-stationary effect of gravitational bending of light caused 
by large clumps of invisible (dark) matter in our Galaxy. Apparent positions of all quasars 
in the sky are shifted from their true positions due to this effect. Since the Sun is moving 
through the Galaxy the angle of the gravitational deflection of light will gradually change 
because of the change in the amount of the matter concentrated near the line of sight towards 
the quasar. The field of the spurious proper motions generated by this gravitational effect 
will look similar to that given by Eq. (34). In other words, the gravitational defiection of 
light can, in principle, interfere with the secular aberration effect. However, according to 
our estimates, the magnitude of this effect is much less than 1 fias and can be neglected. 
Nevertheless, more precise study of the gravitational lensing by clumps of dark matter in 
our Galaxy would be desirable (see, for example, Larchenkova & Doroshenko (1995)). 

The research described in this paper was in part carried out at the Jet Propulsion Lab- 
oratory, California Institute of Technology, under a contract with the National Aeronautics 
and Space Administration. Grant support from the Eppley Foundation for Research and 
Research Council of the University of Missouri-Columbia is greatly appreciated. 
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Fig. 1. — Expected pattern of the proper motion field induced by the secular acceleration 
of the solar barycenter. The length of a 4 //as/yr proper motion is shown in the upper left 
corner. 
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Fig. 2. — Expected pattern of the proper motion field induced by tfie residual rotation of the 
coordinate system about X axis. The rotational rate is made comparable with the amplitude 
of the secular aberration. 



-19- 



80 
60 
40 

D) 
CD 

-D 20 



^^ OH I 11 I t V • 1 1 1 I I I I I I I I f ■ . M l 

^ lllVlvv*, »rrrriiiiiii** ^ 4 J i I I I I 

-= -20 1 I V V V V V *' r r r ! f I \ | | v 1 ^ > ^ ^ ^ . . 7 / i I 

-40 1 V V V V ^ ^ ^ / I I I V \ \ > ^ ^ ^ y / I I 

1 V \) ^ >i P ff 9 « ^ !^ <^ ^ i \ 

-60 I V V ^ ^ ^ / r I \ \ ^ ^ ^ ✓ / i 

-80 ^ ^ t \ ^ ' 



_J_ 



-150 -100 -50 50 100 150 

gal. longitude, deg 



Fig. 3. — Expected pattern of the proper motion field induced by tfie residual rotation of the 
coordinate system about Y axis. The rotational rate is made comparable with the amplitude 
of the secular aberration. 
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Fig. 4. — Expected pattern of the proper motion field induced by tfie residual rotation of the 
coordinate system about Z axis. The rotational rate is made comparable with the amplitude 
of the secular aberration. 
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Table 1. Covariances of the secular acceleration components 
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